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1.0 


INTRODUCTION 

In  this  Final  Technical  Report  the  reaults  of  a  three  yeara* 


research  effort  directed  toward  the  theoretical  study  of  noise  in 
crossed-field  beams  is  presented.  The  major  findings  of  this  study  are: 

(1)  The  only  noise  growth  mechanism  identified  by  the  model 
in  the  gun  region  of  a  crossed-field  beam  is  that  due 
to  the  diocotron  effect. 

(2)  In  order  to  reduce  the  noise  arising  in  the  gun  region  of 
a  crossed-field  beam,  it  is  Important  that  the  entire 
cathode  be  operated  in  the  temperature  limited  regime 
with  focusing  and  other  electrodes  designed  to  give 
smooth,  noncydolding  beams. 

A  more  complete  summary  of  the  objectives  and  major  accomplishments 
of  this  study  is  given  in  the  next  two  sections.  In  the  main  body  of  the 
report,  more  detailed  discussions  of  the  two-  and  three-dimensional  com¬ 
puter  models  are  given  together  with  detailed  discussions  of  the  noise 
simulation  experiments  which  support  the  major  findings  given  above. 


1.1  Research  Objectives 

The  general  research  objectives  as  given  in  the  contract  were: 

The  contractor  shall  furnish  scientific  effort  needed 
to  conduct  the  following  research: 

a.  Carry  out  a  theoretical  analysis  of  growth  of  noise 
in  cross-field  microwave  tubes  based  on  solving  al¬ 
ternately  and  Iteratively  in  steps  two  systems  of 
equations  -  one  system  describing  the  motion  of  each 
electron  and  the  other  Poisson's  equation  depicting 
the  electrostatic  potential  corresponding  to  an  electron 
concentration  -  with  the  noise  Identified  as  averages  of 
parameters  of  the  motion. 


) 


b.  Computer  aspects  of  (a)  are  to  be  performed  employing 
the  most  efficient  computer  software  and  hardware. 


1.1.1  Specific  Research  Tasks  of  the  First  Year's  Research  Effort 

(1)  Conduct  a  preliminary  noise  study  using  the  existing  two- 
dimensional  model  at  cathode. 

(li)  Generalise  model  under  (i)  by  introducing  new  boundary 
conditions  and  variable  model  size. 

(ill)  Investigate  two-dimensional  noise  models  in  parallel 
plane  diodes. 

(iv)  Generalize  to  investigate  three-dimensional  noise  models. 

(v)  Obtain  quantitative  effects  on  noise  due  to  cathode  length, 
magnetic  field  strength,  cathode  current  density,  feedback- 
to-potentlal-minlmum  noise  mechanism,  and  the  diocotron 
noise  mechanism. 

(vi)  Identify  mechanisms  which  on  the  basis  of  the  theoretical 
studies  imply  reduced  noise  levels. 


1.1.2  Specific  Tasks  of  the  Second  Year’s  Research  Effort 

(1)  Implement  the  three-dimensional  noise  model,  paying  special 
attention  to  program  efficiency. 

(ii)  Initiate  preliminary  studies  of  noise  mechanisms  using  the 
three-dimensional  model. 

(ill)  Test  validity  of  present  approximate  techniques  to  study 
arbitrary  electrode  shapes  in  the  two-dimensional  model. 

(iv)  Continue  studies  to  obtain  further  quantitative  data  on 
the  noise  generation  and  noise  growth  mechanisms  in  the 
two-dimensional  model. 
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(v)  Investigate  the  extent  of  computer  generated  noise  in 
the  above  tasks. 

1.1.3  Specific  Tasks  of  the  Third  Year’s  Research  Effort 

(i)  Perform  a  detailed  study  of  noise  mechanisms  using  the 
three-dimensional  model. 

(ii)  Improve  the  program  efficiency  of  the  three-dimensional 
model  developed  under  this  contract  by  coding  portions 
of  the  program  in  assembly  language. 

(ill)  Include  and  study  the  effect  of  grids  in  the  three- 
dimensional  model. 

(iv)  Study  any  noise  reduction  techniques  resulting  from  the 
study  of  noise  mechanisms . 

1.2  Summary  of  Research  Accomplishments 

The  research  accomplishments  associated  with  the  one-  and  two- 
dimensional  models  are  listed  below.  All  of  the  specific  tasks  required 
in  the  first*  second*  and  third  years'  effort  were  completed  with  the 
exception  of  1.1.3  (ii)*  where  time  did  not  permit  coding  of  the  program 
In  assembly  language  to  Improve  program  efficiency*  and  1.1.3  (ill)* 
where  time  constraints  only  made  it  possible  to  study  the  effects  of 
grids  in  one  three-dimensional  gun  geometry.  The  study  of  more  three- 
dimensional  gun  configurations  with  grids  Included  would  certainly 
have  been  desirable. 


It  should  be  noted  that  the  two-  and  three-dimensional  computer 
models  are  not  only  capable  of  giving  information  about  the  noise  generated 
in  a  crossed-field  beam  but  can  also  be  used  to  analyze  and  design  crossed- 
field  guns  in  order  to  obtain  optimum  optics  and  other  beam  characteristics 


Research  Accomplishments  Associated  with  the  Two-Dimensional  Model 

(1)  An  existing  self-consistent  two-dimensional  model  was 
successfully  modified  and  extended  to  permit  the  study  of 
noise  in  two-dimensional  crossed-field  geometries. 

(2)  A  long  Kino  type  crossed-field  gun  was  studied  using  the 
two-dimensional  model.  The  steady  state  characteristic' 
of  the  gun  agreed  closely  with  the  theoretical  values 
predicted  by  the  Kino  gun  theory. 

(3)  Noise  studies  of  the  two-dimensional  Kino  gun  led  to  th 
following  conclusions: 

-  Shot  noise  at  the  cathode  is  the  most  important  con¬ 
tributor  to  gun  noise. 

-  The  noise  is  concentrated  on  the  upper  edge  of  the  beam 
as  it  exits  the  gun  region. 

-  Designing  guns  with  good  beam  optics  appears  to  be  most 
Important  in  reducing  beam  noise. 

-  The  noise  Increases  in  the  gun  region  at  an  exponential 
rate  consistent  with  the  diocotron  growth  mechanism. 


Research  Accomplishments  Associated  with  the  Three-Dimensional  Model 

(1)  A  three-dimensional  conputer  model  was  developed  with  the 
following  important  characteristics: 

-  The  model  is  self-consistent  and  convergent. 

-  The  model  Includes  the  effects  of  variable  magnetic 
fields  in  all  three  dimensions. 


The  model  includes  the  effects  of  space  charge. 


The  model  Is  capable  of  handling  a  variety  of  boundary 
conditions. 


-  The  effect  of  grids  can  be  taken  Into  account  In  the 
model. 

(2)  In  addition  to  verifying  the  conclusions  of  the  two-dimensional 
work,  three-dimensional  noise  studies  of  parallel  plane  diode 
and  Kino  type  guns  have  led  to  the  following  conclusions: 

-  A  small  magnetic  field  perpendicular  to  the  cathode 
does  not  decrease  beam  noise. 

-  The  smoothing  effects  of  the  space  charge  minimum  at 
the  cathode  are  present  and  are  consistent  with  the 
results  obtained  by  Harker. ^ 


1.3  Publications  Resulting  from  the  Work  Performed 

As  a  result  of  the  work  performed  under  this  contract,  a  paper  was 
presented  at  the  I EDM  Meeting  on  December  7,  1977  In  Washington,  D.C. 

A  modified  version  of  this  paper  was  published  in  the  IEEE  Transactions 
on  Electron  Devices.  Both  papers  are  listed  below. 

(1)  E.  K.  Shaw  and  G.  P.  Kooyers,  "Computer  Aided  Design  and 
Analysis  of  Electron  Guns  for  Injected  Beam  Crossed  Field 
Amplifiers,"  presented  at  the  International  Electron  Devices 
Meeting,  Wash.,  D.C.,  1977.  Published  in  the  Technical 
Digest  of  the  IEDM,  p.  440,  1977.  (See  Appendix  B  of  this 
report  for  the  complete  text) 

(2)  E.  K.  Shaw  and  G.  P.  Kooyers,  "Computer  Aided  Design  and 
Analysis  of  Electron  Guns  for  Injected  Beam  Crossed  Field 
Amplifiers,"  IEEE  Transactions  on  Electron  Devices,  Vol. 
ED-26,  No.  7,  p.  1100,  July  1979.  (See  Appendix  C  of  this 
report  for  the  complete  text) 


1.4 


Principal  Investigators 


(1)  Gerald  P.  Kooyers,  President  of  UCA  Systems  Inc. ,  Palo  Alto, 
California. 

(2)  Elden  K.  Shaw,  Senior  Research  Scientist,  UCA  Systems  Inc., 
Palo  Alto,  California*  and  Professor  of  Electrical  Engineering, 
San  Jose  State  University,  San  Jose,  California. 
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HISTORICAL  REVIEW  OF  THE  CROS SED-FIELD  NOISE  PROBLEM 

In  the  first  section  below  the  state-of-the-art  of  crossed-field 


noise  studies  prior  to  1977  is  discussed,  and  the  limitations  and  problems 
associated  with  these  studies  are  pointed  out.  In  the  subsequent  section 
a  comparison  of  the  studies  sponsored  by  the  Air  Force  Office  of  Scientific 
Research  is  given  based  on  1978-1979  annual  reports  and  verbal  conmuni- 
catlon  with  the  other  researchers. 


2.1  The  Problem  of  Noise  in  Crossed-Field  Devices  Prior  to  1977 

Although  a  considerable  amount  of  experimental  and  some  theoretical 
work  was  done  on  noise  in  crossed-field  devices  prior  to  1977,  many  of  the 
fundamental  noise  mechanisms  were  not  well  understood.  Almost  no  quanti¬ 
tative  results  were  available  for  predicting  the  amount  of  noise  generated, 
and  no  definitive  experiments  had  been  made  which  showed  how  to  reduce  the 
noise.  The  justification  for  the  above  conclusions  is  given  in  Appendix  A, 
where  a  review  of  the  noise  problems  in  crossed-field  beams  is  given. 

Ultimately  almost  all  the  noise  in  crossed-field  beams  can  be 
traced  to  fluctuations  in  velocity,  position,  and  emission  time  of  the 
electrons  emitted  off  the  cathode  surface.  A  study  of  noise  thus  must 
begin  at  the  cathode  surface.  If  these  noise  fluctuations  generated  at 
the  cathode  did  not  grow  as  the  electrons  move  downstream,  the  problem  of 
noise  in  crossed-field  beams  would  be  greatly  reduced.  However,  as  dis- 
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cussed  in  Appendix  A,  two  mechanisms  of  noise  growth  have  been  identified, 
and  others  may  exist  when  the  three-dimensional  nature  of  the  beam  is 
considered. ^  In  fact,  one  of  these  noise  mechanisms  (feedback  to  the 
potential  minlmumv  ')  leads  to  a  prediction  of  instability  for  the  "long" 
crossed-field  guns  now  in  use.  Any  attempt  to  model  and  predict  noise 
must  then  Include  a  reasonable  downstream  portion  of  the  beam  and  must 
not  be  limited  to  the  small  signal  regime. 

The  statistical  or  Monte  Carlo  approach  to  the  study  of  noise  in 
crossed-field  guns  was  used  by  Wadhwa  and  Rowe,^  Pollack  and  Whinnery,  ^ 
and  most  recently  by  Lele  and  Rowe.^  Basically  each  of  these  approaches 


is  the  same:  the  beam  la  approximately  almulated  by  roda  of  charge 
(two-dimensional  analyses  only)  which  are  emitted  off  the  cathode  ac¬ 
cording  to  the  random  emission  processes  at  the  cathode  (with,  in  general, 
some  approximations).  The  motion  of  these  rods  of  charge  Is  then  followed 
In  time  with  the  motion  governed  by  the  Lorentz  Force  Law  with  the  electric 
fields  obtained  from  a  solution  of  Poisson's  Equation. 

Unfortunately,  there  are  two  problems  associated  with  each  of 
these  previous  models  which  drastically  limited  the  applicability  and 
usefulness  of  the  results  obtained.  The  first  problem  has  to  do  with 
the  choice  of  beams  studied.  Only  "short"  guns  with  relatively  low  current 
density  were  modeled.  One  of  the  most  Important  noise  mechanisms,  feedback 
to  the  potential  minimum,  thus  was  not  studied  and  high  current  densities 
which  turn  out  to  be  desirable  for  essentially  all  present  CFA's  were  not 
treated  at  all.  The  second  problem  In  each  previous  Monte  Carlo  approach 
has  to  do  with  computer  time,  and  thus  cost.  Because  the  computer  code 
for  the  solution  of  each  case  took  a  relatively  long  time  to  run,  only  a 
very  limited  number  of  cases  was  studied.  The  variation  of  beam  noise  as 
various  parameters  change  Is  not  available  from  these  previous  studies. 

2.2  Comparison  of  AFOSR  Crossed-Fleld  Noise  Studies 

In  Table  I  a  comparison  of  the  four  noise  studies  sponsored  by 
AFOSR  is  presented.  This  table  Is  based  on  1978-1979  reports  ^ 

and  private  communications  from  the  researchers Involved  in  the 
studies.  The  comparison  Is  not  meant  to  be  comprehensive.  Rather,  only 
the  basic  methods  of  approach  and  preliminary  results  are  compared. 


TABLE  I 

COMPARISON  OF  AFOSR  CROSSED-FIELD  NOISE  STUDIES 


Shaw  6 
Kooyers 

UCA 

Systems 


Fontana 
MacGregor 
4  Rowe 
Harris 


Lagranglan  Eulerlan 
(particle)  (fluid) 

Lagranglan n , 
(particle) U‘ 


Harker  & 
Crawford 

Stanford 

Unlv. 

Shkarofsky 

MPB 

Technologies 

Eulerlan 

(fluid) 

Eulerlan 

(fluid) 

Tes 

Yes 

No 

No 

rhree-Dlmenslonal  Model 
-  Completed 


Self  Consistent 
p roach 


Two-Dimensional  Noise 
Calculations 


Verified  Space  Charge 
Smoothing 


Verified  Diocotron 
Growth  Mechanism 


Capable  of  Modeling 
Grids 


Detected  Feedback  to 
Potential  Min.  Instab. 


Capable  of  Modeling 
Entire  Gun  Region 


Useful  as  a  Design 
Tool  for  CFA  Guns 


Three-Dimensional  Noise 
Calculations 


NOTES:  (1)  Fontana  used  a  fluid  approach  In  the  two-dimensional  model  which 
gave  unsatisfactory  results. 

(2)  Both  Harker  and  Shkarofsky  assumed  a  parabolic  d.c.  potential  at 
the  cathode. 

(3)  Harker  solved  for  current  fluctuations  only. 

(4)  Shkarofsky  solved  for  velocity  fluctuations  only. 
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No 

No 

No 

i 

No 
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THE  TWO-DIMENSIONAL  NOISE 


In  the  first  year  s  study  a  two-dimensional  model  which  was 
originally  used  as  a  design  tool  for  crossed-field  guns  was  refined 
and  modified  for  the  study  of  noise  in  crossed-field  guns.  The  details 
of  this  two-dimensional  siodel  are  given  in  Appendices  B  and  C  of  this 
report.  In  Fig.  .3.1  a  flow. diagram  is  presented  of  the  overall  approach 
used  in  both  the  one-  and  two-dimensional  studies. 

The  method  of  simulating  the  noisy  beam  at  the  cathode  and  the 
method  of  calculating  noise  from  the  beam  characteristics  are  discussed  in 
detail  below  since  these  topics  are  not  discussed  in  detail  in  the  Appendices 


.1  Simulation  of  the  Noisy  Beam  at  the  Cathode 


Since,  as  previously  stated,  the  ultimate  source  of  nearly  all 
the  noise  is  due  to  random  processes  of  electron  emission  at  the  cathode 
surface,  it  is  very  Important  that  the  emission  properties  of  the  cathode 
are  modeled  as  accurately  as  possible.  The  cathode  emission  is  completely 
specified  by  the  following: 

(1)  The  number  of  electrons  emitted  in  a  time  interval  At 

(2)  The  emission  velocities  of  the  electrons 

(3)  The  time  of  emission  of  the  electrons. 

The  method  for  calculating  each  of  the  above  three  quantities  is  summarized 
below. 


3.1.1  Emission  Number 

It  can  be  shown that  the  probability  that  exactly  S  electrons 
are  emitted  from  a  given  cathode  spot  of  length  Ay  during  the  time  in¬ 
terval  At  is  given  by 


I 


Figure  3-1.  Nolee  Simulation  and  Evaluation  Flow  Diagram 


where  X  is  the  average  number  of  electrons  emitted  in  each  At.  This 
is  recognized  to  be  the  Poisson  Distribution. 

The  probability  that  S  or  fewer  electrons  are  emitted  in  the 
time  interval  At  is  given  by  the  cumulative  distribution  function 

S 

F(S)  -  l  f(S)  (3-2) 

S-0 

and 

0  <  P(S)  *  1  (3-3) 

To  determine  the  number  of  electrons  emitted  in  a  time  interval*  At* 
a  uniformly  distributed  random  number*  R»  between  0  and  1  is  generated. 
The  number  of  electrons  emitted  is  S  if 

F(S-l)  <  R  «  F(S)  (3-4) 

When  R  .<  f(0)  no  electrons  are  emitted. 

For  computational  purposes*  the  following  procedure  is  then 

followed: 

(1)  The  maximum  number*  N  *  of  electrons  allowed  to  be 

1& 

emitted  from  the  area  AyAz  and  in  any  At  is  selected. 

(2)  f(S)  is  then  calculated  for  each  value  of  S  between 
0  and  M  . 

B 

(3)  F(S)  and  F(S-l)  are  then  calculated  for  each  S  between 
0  and  N  using  the  calculated  density  functions  f(S) 
and  are  then  stored  in  a  table  containing  S*  F(S),  and 
F(S-l) • 

(4)  Random  number  R  is  then  generated. 


(5)  Knowing  R,  the  pair  of  values  P(S-l)  and  7(S)  that 
bracket  E  are  found  fro*  the  table  and  the  corresponding 
value  of  S  la  then  known. 

Initially  a  value  of  20  was  chosen  for  This  value  was  used  by 

O' Flynn  and  Kooyers^^  for  siallar  problems. 


3.1.2  Emission  Velocities 

A  half-Maxwelllan  distribution  Is  assumed  for  the  x  component 
of  emission  velocity.  The  x  velocity  Is  then  given  by ^ 


(“In  Rx)* 


(3-5) 


where  R  Is  a  random  number  uniformly  distributed  between  0  and  1, 

T  Is  the  cathode  temperature,  k  is  Boltzmann's  constant,  and  m  is 
the  electron  mass. 


A  full-Maxwelllan  distribution  Is  assumed  for  electron  velocities 
In  the  y  direction  (for  the  two-dimensional  case)  and  the  y  and  z 
directions  (for  the  three-dimensional  case).  The  relation  giving  these 
velocities  Is 


e 

y 


erf*1(Ry) 


(3-6) 


and 


z  • 


2kT 


.rf-2«,> 


(3-7) 


where  R 

y 

and  1.  In 
formula  is 


and  R  are  random  numbers  uniformly  distributed  between  0 
s 

order  to  simulate  these  velocities  efficiently,  an  approximate 
used  for  the  error  functions. 
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3.1.3  &!<«•< on  Times 

In  order  to  avoid  calculating  the  positions  of  the  electrons 
above  the  cathode  for  random  emission  times,  we  assume  that  all  S  of 
the  electrons  in  a  given  Ay  are  emitted  at  the  same  time  but  with  a 
random  distribution  over  the  by.  The  emission  points  are  then  given  by 

y  ■  ±  Ry  (3-8) 


t 


m  ±  V  ^1 
*0  *  2 


(three-dimensional  case) 


(3-9) 


where  R  is  a  random  number  between  0  and  1  and  the  sign  of  R  is 

y  y 

also  selected  randomly.  Thus,  to  an  observer  looking  at  one  point  on  the 

cathode  surface,  the  electrons  will  appear  to  come  off  that  point  at 

random  times.  Provided  the  time  interval  bt  is  chosen  to  be  relatively 

small,  the  degree  of  approximation  introduced  should  be  minor. 


3.2  Noise  Analysis  Procedures 

From  the  computer  model  we  have  available  the  positions  and 
velocities  of  each  electron  at  each  discrete  time  t  •  NfcAt  where  Nfc 
is  the  iteration  number.  In  order  to  be  useful,  this  information  must 
be  related  to  the  noise.  There  are  several  different  noise-related 
parameters  which  are  calculated  in  this  study. 

In  general  we  would  like  to  be  able  to  calculate  each  of  these 
parameters  at  any  arbitrary  point  (x,y,c)  within  the  beam.  However,  the 
nature  of  the  way  we  solve  the  problem  using  a  finite  number  of  electrons 
limits  us  to  finding  the  fluctuations  in  some  small  but  finite  region. 

In  order  to  obtain  enough  "electrons"  in  each  region  for  reasonably  good 
statistics,  the  gun  region  is  divided  into  512  noise  regions  which  are 
3x3  mesh  units  in  size.  These  regions  are  sketched  in  Figure  3-2. 

The  method  used  for  calculating  each  of  the  parameters  of  interest  is 
discussed  below. 
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Figure  3-2.  Noise  Regions  in  the  CFA  Gun  Region. 
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3.2.1  A.C.  Velocity  Components 

In  order  to  obtain  the  a.c.  velocity  components  we  first  must 
find  the  average  velocity.  The  a.c.  velocity  Is  then  equal  to  the  actual 
velocity  minus  the  average  velocity.  The  velocity  averaged  over  one  noise 
region  at  a  given  tine  t  is  given  by 


N. 


i  V 

l  v  (t.i) 
l.J  £-1  xi,J 


(3-10) 


where  v  (t,£)  Is  the  x  velocity  of  the  £th  electron  in  the 

l.J 

1, J-th  noise  region.  Since  the  equations  for  the  other  component 
are  of  the  sane  form,  ve  carry  the  analysis  only  for  the  component. 

Over  the  tine  Interval  T  •  N^At  the  time  average  velocity  Is  given  by 


K 


N 


The  a.c.  velocity  component  can  now  be  calculated  at  time  t  from 
V  (t)  -<V 

l.J  i.J  l.J.k 


(3-11) 


(3-12) 


where  v  (t)  Is  the  a.c.  component  of  velocity  in  the  l,J-th  noise 

Xi.J 

region  at  time  t. 


3.2.2  Variance  of  the  Velocity  Components 

The  variance  of  the  x-.eloclty  at  time  t  can  be  found  from 


l.J 


(t)  -  <vj  <t)>  -  <v  (t)>- 


(3-13) 


l.J 


l.J 


I; 


where  o2  (t)  Is  the  variance  of  the  x-component  of  velocity  In  the 

*!»  J 

i,J-th  noise  region  at  time  t  and  <v2  (t)>  Is  given  by 

*i.J 


» 
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(3-14) 


1  J 

(t)^--  — —  {  vj  (t,£) 

*1,J  W1,J  £-1  Xl,  j 


The  variance  of  the  other  component  of  velocity  can  be  found 
In  a  alailar  Banner.  The  variance  and  the  a.c.  components  of  velocity 
are  useful  In  that  they  give  the  relative  magnitude  of  the  fluctuations. 
Thus  the  "growth"  of  these  fluctuations  can  be  studied  by  calculating 
the  a.c.  velocity  components  and  variances  at  various  downstream  posi¬ 
tions  In  the  besa.  The  a.c.  cosponents  of  velocity  and  their  variance 
do  not  give  any  Information  about  the  frequency  spectrum  of  these 
components. 


3.2.3  A.C.  Current  Components 


The  charge  density  In  the  i, j-th  noise  region  at  time  t  Is 


given  by 


±1  "o  1 

AxAy 


(3-15) 


where  Is  the  scaling  factor  for  the  electrons.  The  time  averaged 

space  charge  density  in  the  l,J-th  noise  region  Is  given  by 


"  AxAyAz°  Nj. 


(3-16) 


Mow  the  current  density  at  time  t  in  the  1, j-th  rectangle  Is  given  by 


J  (t)  -  p.  ,(t)  x  .(t) 

*i,j  1.3  1.3 


(3-17) 


and  the  time  average  current  density  Is  given  by 


<J*  >  "<pi,j>  *<Vx  > 

*i.j  1,3  *1.3 


(3-18) 
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The  a.c.  component  of  the  x- component  of  current  density  Is  thus  given 
hy 


Jx  (t)  -  (t)  -<J  > 

*i»J  i.J 


(3-19) 


Again,  although  this  parameter  does  not  tell  us  anything  about  the  fre¬ 
quency  components  of  the  fluctuations.  It  Is  very  useful  in  determining 
the  relative  growth  rate  of  current  fluctuations  In  the  beam. 


3.2.4  Correlation  Functions  and  Power  Spectra 

Finally,  In  order  to  completely  determine  the  noise  properties 
of  the  beam.  It  Is  necessary  to  calculate  the  relative  magnitudes  of 
the  various  frequency  components  of  the  noise  at  a  particular  cross- 
section  of  the  beam.  The  two  quantities  of  interest  In  beginning  this 
calculation  are  the  a.c.  kinetic  potential  and  the  a.c.  current  fluctuation. 
These  two  quantities  may  be  defined  for  a  two-dimensional  beam  as^ 


V(t)  - 


-jiQ  •  v(t) 


(3-20) 


J(t)  - 


3(t)  •  X 


(3-21) 


where  V(t)  Is  the  a.c.  kinetic  potential  at  time  t 
Is  the  d.c.  velocity  of  the  beam 
v(t)  Is  the  a.c.  component  of  the  velocity  of  the  beam 
n  Is  the  charge-to-ma8s  ratio  of  the  electron 

J(t)  Is  the  a.c.  current  density  fluctuation  at  time  t 

3(t)  is  the  a.c.  current  density  in  the  beam 

Is  the  d.c.  current  density  in  the  beam. 
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At  this  point  the  final  problem  in  the  spectral  analysis  is 
in  relating  a  spectrum  consisting  of  finite  discrete  data  of  a  function 
V(t)  or  J(t)  (each  value  obtained  as  an  average  over  a  small  Interval 
At)  to  the  true  spectrum  of  the  continuous  time  function.  Several  authors 
have  treated  this  subject  in  some  detail.  Results  show  that  it  is  possible 
to  calculate  the  Haus  noise  parameters  8(f),  ♦  ( f),  and  8(f)  from  the 

discrete  data  where 

0(u>)  is  the  self  power  density  spectrum  of  the  kinetic  noise 
voltage, 

8(u>)  i»  the  self  power  density  spectrum  of  the  noise  current 
modulation,  and 

0(u)  is  the  cross-power-density  spectrum  between  the  kinetic 
voltage  and  the  current  modulation. 

For  a  crossed-field  beam  an  exact  relation  between  the  noise 
parameters  and  the  noise  figure  has  not  been  derived.  However,  even 
without  an  exact  expression,  a  study  of  the  relative  magnitudes  of  the 
various  noise  parameters  (0,  ♦,  and  8)  as  a  function  of  position  and 
other  gun  and  beam  parameters  is  very  valuable  in  gaining  an  understanding 
of  the  noise  mechanisms  in  crossed-field  electron  devices. 


4.0 


TWO-DIMENSIONAL  NOISE  SIMULATION  EXPERIMENTS 


Using  the  two-dimensional  model,  several  noise  simulation 
experiments  have  been  carried  out  on  a  long  Kino  gun.  The  geometry 
and  parameters  of  the  long  Kino  gun  are  given  In  Appendix  C.  In  parti¬ 
cular,  the  noise  growth  rate,  the  variation  of  noise  throughout  the  gun 
region,  and  the  effect  of  changes  in  the  beam  injection  conditions  have 
been  studied. 


4.1  Noise  Growth  in  the  Gun  Region 

In  Figure  4-1  the  noise  power  spectral  density  due  to  fluctuations 
in  the  kinetic  potential  is  plotted  from  the  end  of  the  cathode  to  where 
the  beam  exits  the  gun  region.  The  noise  growth  appears  to  be  approxi¬ 
mately  linear  when  plotted  on  semi-log  paper,  and  thus  the  actual  noise 
growth  is  exponential.  For  the  case  of  Figure  4-1  a  straight  line 
approximation  gives  a  noise  growth  of  approximately  170  db/in  which  is 
a  rather  substantial  growth  rate.  The  diocotron  theory  also  predicts 
an  exponential  growth  of  the  noise. 


4.2  Noise  Variation  Across  the  Beam 

Shown  in  Figure  4-2  and  Figure  4-3  is  the  variation  in  noise 
from  the  bottom  to  the  top  of  the  beam  as  the  beam  exits  the  gun  region. 

In  all  of  these  caseB,  noise  is  greater  on  the  top  edge  of  the 
beam  then  on  the  bottom  edge.  This  may  be  due  to  the  fact  that  the 
electrons  on  the  top  edge  of  the  beam  have  traveled  a  greater  distance 
and  thus  have  experienced  more  diocotron  nol6e  growth  than  those  electrons 
on  the  bottom  edge  of  the  beam. 

Unfortunately,  as  the  beam  enters  the  interaction  region  of  a 
CFA,  the  noisiest  part  of  the  beam  will  be  nearest  the  slow  wave  circuit 
and  thus  will  be  more  strongly  coupled  to  the  slow  wave  circuit. 
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Variation  of  Noise  With  Frequen 


Shown  in  Figure  4-4  and  Figure  4-5  is  the  variation  in  bean 
noise  as  a  function  of  frequency  for  two  different  positions  in  the  bean. 
For  Figure  4-4  the  noise  is  sanpled  near  the  cathode;  while  in  Figure 
4-5  the  noise  is  plotted  for  a  noise  region  near  the  gun  exit. 

At  the  gun  exit*  the  noise  as  neasured  by  the  self  power  spectral 
density,  p,  of  the  kinetic  potential  is  fairly  constant  with  frequency, 
while  the  self  power  spectral  density  of  the  current  density,  p,  has 
large  peaks  and  valleys  which  roughly  correspond  to  nultiples  of  the 
cyclotron  frequency. 

Near  the  cathode  both  p  and  p  vary  with  frequency  and  have 
roughly  the  sane  valleys  and  peaks  corresponding  to  nultiples  of  the 
cyclotron  frequency.  The  variation  of  noise  with  frequency  near  the 
cathode  would  be  unexpected  if  it  were  not  for  the  fact  that  many  of 
the  electrons  near  the  cathode  have  been  in  the  gun  for  a  relatively 
long  time  due  to  two  effects: 

(1)  Electrons  returning  to  the  gun  on  a  cycloidal  path  from 
further  upstream 

(2)  Electrons  which  migrate  downstream  while  maintaining  a 
position  very  close  to  the  cathode. 

In  both  cases  noise  growth  mechanisms  have  sufficient  time  to  affect  the 
noise  characteristics  of  these  electrons. 


Variation  of  Noise  with  Cathode  Entrance  Conditions 


Relative  Frequency  (F/Fc) 


Fig.  4-4.  Self  Power  Spectral  densities  for  Current  and  Kinetic 
Potential  Fluctuations  Near  the  Cathode. 


Relative 

Self 

Power  .8 

Spectral 

Density 

(XIO'*) 


.6 


(1)  Random  positions  and  random  numbers*  Injected  but  with 
four  different  injection  velocities: 

(a)  Zero  velocity  (Temperature  ■  0) 

(b)  Uniform  velocity 

(c)  Maxwellian  velodtlea  (Temperature  ■  T) 

(d)  Velocities  consistent  with  the  electric  fields  at 
the  cathode. 

*  But  with  constant  average  current  emitted  from  the 

cathode  surface  and  with  the  cathode  space  charge  limited. 

(2)  Random  positions,  random  numbers  injected,  and  Maxwellian 
velocities  but  with  variable  amounts  of  charge  per  "electron" 
from  the  cathode  surface. 

Noise  outputs  for  the  various  input  conditions  of  cases  (1)  above 
indicate  that  the  noise  is  not  appreciably  affected  by  the  electron  injection 
velocities.  A  typical  plot  of  noise  output  for  simulated  cathode  tempera¬ 
tures  of  0  and  1323°K  is  given  in  Figure  4-1. 

In  all  the  cases  (1)  above,  the  cathode  is  running  In  a  space- 
charge  -limited  condition.  If  the  numbers  of  electrons  emitted  from  the 
cathode  are  reduced  to  the  point  where  the  cathode  becomes  temperature 
limited,  the  optics  of  the  gun  deteriorate  and  the  noise  in  the  beam 
increases  dramatically. 

When  the  total  charge  carried  by  one  "electron"  is  increased, 
the  shot  noise  is  increased  in  direct  proportion.  This  corresponds  to 
the  input  condition  of  cases  (2)  above.  Noise  output  from  the  computer 
simulation  as  measured  by  the  velocity  variance  also  goes  up  in  direct 
proportion  to  the  amount  of  charge  per  "electron."  These  results  indicate 
that  shot  noise  is  the  main  contributor  to  beam  noise  in  this  model. 


5.0 


THE  THREE-DIMENSIONAL  NOISE  MODEL 


The  three-dimensional  nolee  model  eolvee  for  the  notion  of  charges 
in  e  three-dimensional  rectangular  geometry  Including  the  effects  of  static 
electric  and  nagnetlc  fields  and  opace  charge.  A  computer  program  written 
in  FORTRAN  IV  carries  out  the  steps  given  in  the  flow  chart  of  Figure  5-1. 

A  sample  of  the  output  generated  by  this  program  is  given  in 
Appendix  F.  Information  is  printed  out  giving  cathode  current  density, 
distribution  of  charges  in  the  beam,  selected  charge  trajectory  information, 
etc.  for  each  time  step. 

Belov,  details  are  given  on  the  potential,  trajectory,  and  noise 
calculations. 

The  three-dimensional  program  consists  of  the  major  program 
modules  shown  in  Figure  5-2.  The  function  of  each  program  module  is  given 
below: 


GUN 3D  - 

Main  calling  program  which  reads  in  gun  data, 

initializes  program  variables,  and  calls  other 

program  modules. 

T  RAJ  3D  - 

Solves  the  three-dimensional  trajectory  equations 

including  variable  magnetic  fields  in  the  x,  y, 

and  z  directions. 

EMIT3D  - 

Calculates  particle  emission  numbers,  positions, 

and  velocities  as  they  are  emitted  from  the  cathode 

P0SS3D  - 

Solves  Poisson's  Equation  in  three  dimensions  using 

subroutines  FOUR67  and  PUSCRT. 

N01S3D  - 

Calculates  current  and  velocity  fluctuations  in  the 

beam  and  uses  these  fluctuations  to  calculate  the 
space  charge  smoothing  factor  and  the  "temperature" 
of  the  beam  at  various  locations  within  the  beam. 
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Figure  5-1.  Flow  Chart  Showing  Major  Steps  in  the 
Three-Dimensional  Program. 


Solving  Poisson's  Equation  In  Three  Dimensions 


In  order  to  obtain  the  potentials  in  the  three-dimensional  region, 
it  is  necessary  to  solve  the  Poisson  Equation 


524>  32<t 

3*2  +  awz 


(5-1) 


In  the  region  sketched  in  Figure  5-3.  The  boundary  conditions  shown  in 
Figure  5-3  are  chosen  to  be  consistent  with  three-dimensional  crossed- 
fleld  gun  geometries.  Boundary  conditions  on  the  z-boundaries  allow  for 
potentials  with  the  derivative  equal  to  zero,  periodic  potentials,  or 
constant  potentials. 

The  finite  difference  techniques  used  In  the  solution  are  given 
in  detail  in  Appendix  D  for  the  case  of  fixed  potentials  on  the  z-boundaries 
and  follow  a  method  suggested  by  Svartrauber  and  Sweet.  Figure  5-4 

shows  the  mesh  geometry  used  in  the  finite  difference  solution.  A  very 
brief  review  of  the  technique  is  given  below. 

(1)  The  finite  Fourier  expansion  of  the  five  point  difference 
approximation  to  Poisson's  Equation  is  taken  in  the 
z-direction  using  a  program  which  performs  the  Fourier 
transform^^  (FOUR67) .  This  leads  to  n-1  two-dimensional 
Helmholtz  equations 


V2i  +  K2$  -  p 


(5-2) 


where  4  and  p  are  the  transformed  potentials  and  space 
charge  and  K  is  a  constant  which  depends  on  n.  The 
space  charge  assigned  to  each  mesh  point  is  found  by  adding 
together  the  total  number  of  particles  in  a  cell  (hx  i  x  h^) 
centered  on  the  mesh  point.  This  is  commonly  called  the 
particle  in  cell  (PIC)  approach. 

(2)  Equation  (5-2)  is  then  solved  for  the  transformed  potentials 
using  an  existing  program  (PWSCRT) . 


31 


(x,y,c)  -  0 


zero  derivative 

periodic 

constant 


Bz 

♦(x,y,c)  -  *(x,y,0) 


|f(x,y,0)  -  0  ■  0 


$(x,y,0)  »  $(x,y,c) 

♦(x,y,0)  -  *c 

Figure  5-3.  Three-Dimensional  Geometry  with. 
Boundary  Conditions. 


n-di visions 


Figure  5-4.  Three-Dimensional  Geometry  with 
Modified  Boundary  Conditions. 


(3)  Given  the  transformed  potentials  a  program  (FOUR67) 

is  used  to  obtain  the  actual  potentials  by  performing  the 
Inverse  Fourier  Transform. 


5.2  Equations  of  Motion 

A  detailed  discussion  of  the  development  of  the  three-dimensional 
finite  difference  equations  of  motion  is  given  in  Appendix  E.  There  are 
several  facts  which  should  be  noted  concerning  these  equations: 

(1)  The  finite  difference  equations  of  motion  are  accurate 
to  second  order  and  are  solved  in  a  right-handed  coordi¬ 
nate  system. 

(2)  The  electric  field  is  assumed  to  be  a  slowly  varying 
function  o.f  position  (i.e.,  does  not  change  much  in  one 
time  step) . 

(3)  The  magnetic  field  is  assumed  to  consist  of  x,  y,  and  z 
components  which  may  be  slowly  varying  functions  of 
position. 

(4)  The  equations  are  solved  in  such  a  way  as  to  be  exact 
for  any  size  time  step  provided  the  electric  and  magnetic 
fields  are  constant. 


5.3  Noise  Analysis  Procedures 

A  separate  noise  analysis  program  (N01S3D)  has  been  written 
to  perform  noise  analysis  on  the  three-dimensional  beam.  This  program 
calculates  two  parameters  which  are  closely  related  to  beam  noise:  the 
space  charge  smoothing  parameter  (current  fluctuations)  and  the  beam 
temperature  (velocity  fluctuations).  Details  of  the  techniques  used  to 
make  these  calculations  are  given  below. 


5.3.1  Current  Fluctuations 

In  order  to  determine  the  growth  of  current  fluctuations  in 
the  beam,  the  current  fluctuations  at  the  cathode  are  calculated  using 
computer  generated  values  of  Instantaneous  cathode  current,  Ig,  and 
a  running  average  of  the  cathode  current,  <Ig> .  At  each  Iteration,  the 
current  fluctuation  at  the  cathode  is  calculated  from 

\  ’  h  -  <V  (5*3> 


Similarly  the  fluctuations  in  the  beam  current  at  various  planes 
in  the  gun  region  are  calculated  from 

i,  -  h  -  'V  (5-‘> 


at  each  Iteration. 


A  space  charge  smoothing  factor,  T,  is  then  calculated  by 
taking  the  r.m.s.  of  the  ratio  of  the  current  fluctuations  for  the  beam 
at  a  given  plane  within  the  gun  to  the  cathode  current  fluctuations. 

The  space  charge  smoothing  factor  is  then  defined  to  be 


(5-5) 


5.3.2  Velocity  Fluctuations 

If  the  beam  can  be  considered  to  be  in  a  state  of  equilibrium 
with  a  Maxwellian  distribution  of  velocities,  then  the  temperature  of 
the  beam  may  be  calculated  from 

T  ■  ?r  <(v  -  <v>)^>  (5-6) 


where  m  is  the  electron  mass,  k  is  Boltzman's  constant,  and  v  is 
the  instantaneous  velocity.  In  terms  of  the  components  of  velocity. 


the  above  equation  may  be  written  aa 


T  “  —  (<v  2>  +  <v  2>  +  <v  2>  -  <v  >2  -  <v  >2  -  <v  >2  ) 

3k  x  y  z  x  y  z 

(5-7) 

where  the  bracketa,  <  >,  represent  spatial  averages  in  the  region  under 
analysis.  Running  averages  (time  averages)  of  the  x,y(z  velocities  in 
each  region  under  analysla  are  calculated  at  each  time  step  and  used  to 
determine  the  instantaneous  velocity  fluctuations. 

The  temperature,  T,  is  calculated  in  rectangular  regions  one 
mesh  distance  thick  in  the  xy-plane,  the  xz-plane,  and  the  yz-plane  as 
sketched  in  Figure  5-5.  It  is  thus  possible  to  compare  the  beam  tempera¬ 
ture  at  various  x,  y,  and  z  planes. 


5.4  Inclusion  of  Electrodes  (Grids)  in  the  Model 

Two  approaches  to  the  Inclusion  of  the  grids  in  the  three-dimensional 
geometry  are  discussed  below: 


5.4.1  Approximate  Technique  Using  D.C.  Potentials 

In  this  approach  the  principle  of  superposition  is  used  to  Include 
Internal  electrodes  in  an  approximate  way  as  follows: 

(1)  Using  the  exact  geometry  of  the  gun  including  grids  and 
other  Internal  electrodes,  the  potential  is  solved  in  the 
three-dimensional  geometry  neglecting  space  charge.  In 
the  discussion  below  these  potentials  are  called  the  d.c. 
potentials . 
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Figure  5-5.  Noise  Regions  Used  in  Noise 
Calculations. 


(2)  Poisson  s  Equstlon  Is  solved  for  the  potentials  In  a 
rectangular  "box"  (with  zero  potential  on  the  electrodes) 
which  approximates  the  exact  gun  geometry  as  closely  as 
possible  and  the  cathode  exactly. 

(3)  The  d.c.  potentials  are  added  to  the  space  charge  potentials 
to  obtain  the  potentials  used  In  the  trajectory  equations. 

(4)  Steps  (2)  and  (3)  are  repeated  for  each  Iteration.  Note 
that  the  d.c.  potentials  need  only  be  calculated  once  and 
then  saved. 

A  three-dimensional  relaxation  program  for  solving  Laplace's 
Equation  with  arbitrary  internal  electrodes  and  boundary  conditions  has 
been  written.  The  program  has  been  tested  on  a  three-dimensional  crossed- 
field  diode  with  grids  in  a  16  x  16  x  16  mesh  convergence  point  system. 
After  50  Iterations  (.08  minutes  CPU  time)  the  changes  in  potential  are 
less  than  1Z  and  the  program  appears  to  have  converged. 

This  approach  is  approximate  In  that  the  imaging  of  the  electrons 
In  the  grid  are  neglected.  Since  the  grids  are  normally  run  at  a  negative 
potential,  this  error  Is  expected  to  be  relatively  small. 


5.4.2  Exact  Accroach  Usine  a  Relaxation  Procedure  on  the  Space  Charge 


(18) 

Hockney v  '  has  suggested  several  methods  of  including  the  effect 
of  electrodes  In  the  Interior  of  the  mesh  region  used  for  the  solution  of 
Poisson's  Equation.  When  the  number  of  electrode  points  is  large,  he 
suggests  that  the  space  charge  on  individual  mesh  points  be  adjusted 
iteratively  in  response  to  the  local  error  in  potential  on  that  point. 

This  approach  has  been  progressed  as  a  test  case  in  the  Poisson  solving 
routines.  The  procedure  used  is  outlined  below: 
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(1)  Storage  in  an  arr ay  is  provided  for  the  coordinates,  potential, 
and  space  charge  of  each  internal  electrode  point. 

(2)  An  initial  value  of  charge  is  assigned  to  each  internal 

electrode  point  assuming  the  electrode  is  an  isolated  sphere 

of  dimension  h  X  h  X  h  • 
x  y  z 

(3)  After  each  Iteration  of  Poisson's  Equation  a  new  value  of 
charge  is  assigned  to  each  electrode  point  according  to  the 
following  relaxation  formula 

Q“”  ■  q°id  [l + *«('  pL\  ° )_ 

where  R^  is  a  relaxation  factor  &  1,  Vq  Is  the  electrode 
potential  desired,  and  V  is  the  maximum  of  VQld  and  VQ. 

The  above  technique  has  been  tried  on  a  parallel  plane  geometry 
-test  case  with  one  interior  electrode  point  at  -100  volts.  The  procedure 
was  found  to  converge  with  a  relaxation  factor  of  .25.  Time  did  not  permit 
the  Inclusion  of  this  technique  in  the  study  of  a  gridded  gun.  Obviously, 
this  technique  will  result  in  an  Increase  in  run-time  for  the  computer  program. 
However,  the  increase  in  run-time  may  be  as  small  as  a  factor  of  two  since 
the  technique  can  be  applied  without  a  time  penalty  during  the  iterations 
necessary  for  the  solution  to  reach  a  steady  state.  Once  the  steady  state 
is  reached  the  correction  to  the  charges  necessary  to  maintain  the  correct 
Internal  electrode  potential  should  be  small  and  thus  easily  accomplished 
with  a  few  iterations. 
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6.0 


THREE-DIMENSIONAL  NOISE  SIMULATION  STUDIES 


Using  the  three-dimensional  model,  several  test  cases  were  run 
to  determine  the  correctness  of  the  model  and  obtain  preliminary  noise 
information. 

In  testing  the  model  for  accuracy  of  potential  and  trajectory 
calculations,  a  16  x  16  z  16  mesh  was  used  in  the  parallel  plane  geometry 
of  Fig.  6-1.  For  the  noise  calculations,  a  32  x  48  x  16  mesh  was  U8ed. 

For  the  larger  mesh  system  with  16,000  charges,  each  case  required  approxi¬ 
mately  12  minutes  of  computer  time  on  the  IBM  370/168.  The  results  ob¬ 
tained  using  the  larger  mesh  are  discussed  below. 


6.1  The  Long  Kino  Gun  Test  Case 

As  a  preliminary  test  case,  the  Long  Kino  Gun  sketched  in  Figure 
6-1  was  approximately  modeled  in  three  dimensions  using  a  32  x  48  x  16  mesh 
system.  The  z-boundarles  are  taken  to  be  at  zero  potential  in  order  to 


confine  the  beam.  The  potential  on  the 
as  zero  potential  with  the  potential  on 
from  the  Kino  gun  theory.  This  gun  has 

Total  Beam  Current 
Cathode  Length 
Cathode  Width 
Magnetic  Field 
Current  Density 
Accelerator  Potential 
Normalized  Cathode  Length 


lower  boundary  (cathode)  is  taken 
the  anode  fixed  by  the  potentials 
the  following  design  parameters: 

3.34  amps 

.150  Inches 

.247  inches 

2500  gauss 
2 

14  amp /cm 
5000  volts 
116.4  Kino  units 


Figure  6-2  shows  the  convergence  characteristics  of  the  computer 
simulation.  In  this  figure  the  cathode  current  reaches  a  "steady  state"  in 
about  10  iterations  and  then  fluctuates  about  this  "steady  state"  for  time 
steps  10  through  75.  The  computer  calculated  average  exiting  beam  current 
is  3.62  amps  which  is  near  the  design  value. 
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6.1.X  Current  Fluctuation  Calculations 

Shown  In  Figure  6-3  Is  the  space  charge  smoothing  factor,  I*»  as 
a  function  of  Iteration  number  for  the  Long  Kino  Gun  test  case  with  the 
cathode  operating  space  charge  and  temperature  limited.  In  the  space 
charge  limited  case,  the  current  fluctuations  are  smaller  at  the  beam 
exit  than  at  the  cathode  surface.  This  Indicates  space  charge  "smoothing" 
by  the  potential  minimum  and  is  consistent  with  the  results  obtained  by 
Harker. ^  In  the  temperature  limited  case  the  current  fluctuations  at 
beam  exit  are  larger  than  at  the  cathode  surface.  In  an  ideal  case  of  no 
noise  growth  mechanisms,  the  value  of  T  for  the  temperature  limited  case 
would  be  unity.  The  value  of  T  •  1.75  is  due  to  diocotron  noise  growth 
in  the  beam  as  it  travels  from  the  cathode  to  the  beam  exit  plane.  Space 
charge  smoothing  is  shown  even  more  clearly  in  Figure  6-4  where  the  space 
charge  smoothing  factor,  r2,  Is  plotted  as  a  function  of  the  relative 
current  emitted  off  the  cathode.  The  value  of  T2  is  calculated  at  a 
position  above  the  cathode  where  T2  is  minimum.  As  expected,  the  space 
charge  smoothing  factor  is  approximately  unity  in  the  temperature  limited 
regime  and  drops  dramatically  as  the  cathode  is  operated  space  charge 
limited.  Two  test  cases  were  run  with  a  small  x-directed  magnetic  field. 
The  value  of  T  after  75  iterations  are  given  below: 

Case  I  -  Mo  magnetic  field  in  x-dlrection 

B  •  .25  w/m^ 

z 

r  •  .160 

Case  II  -  B  ■  51  B 

K  2  * 

B^  •  .25  w/m 

T  -  .2005 

Case  III  -  B  -  10Z  Bt 

x  2  * 

B_  ”  .25  w/m 

z 

r  -  .1883  . 

In  both  Cases  I  and  III  the  x-dlrected  magnetic  field  Increases 
the  exiting  current  fluctuations. 


43 


Space 

Charge 

Smoothing 

Factor 

r2 


6.1.2  Bean  Temperature  Calculations 

Shown  In  Figures  6-5 ,  6-6*  and  6-7  is  the  temperature  of  the  beam 
as  a  function  of  noise  region  in  the  x,  y,  and  z  planes.  Figure  6-5  shows 
how  the  bean  temperature  Increases  toward  the  beam  exit.  This  growth 
appears  to  be  nearly  exponential  over  the  cathode  surface,  drops  an  order 
of  magnitude  at  the  end  of  the  cathode,  and  then  rises  again  at  the  same 
exponential  rate.  This  exponential  growth  seems  to  be  due  to  dlocotron 
noise  growth. 

Temperature  variations  as  a  function  of  position  above  the  cathode 
are  shown  in  Figure  6-6.  As  found  in  the  two-dimensional  studies,  the 
top  of  the  beam  seems  to  be  noisier  than  the  bottom.  This  could  be  due 
to  dlocotron  growth  variations  since  the  electrons  on  the  beam  top 
travel  further  than  those  on  the  bottom. 

In  Figure  6-7  the  variation  of  beam  temperature  Is  given  In  the 
direction  of  the  magnetic  field  for  the  case  of  no  magnetic  field  in  the 
x-direction  and  a  small  magnetic  field  in  the  x-direction.  Note  that  the 
magnetic  field  increases  the  temperature  of  the  beam. 


6.2  Three-Dimensional  Studies  With  Grids 

As  a  preliminary  test  case  a  three-dimension  parallel  plane  diode 
was  studied  using  a  16  x  16  x  16  mesh  system.  The  diode  had  a  single  grid 
at  -100  volts  running  the  full  length  of  the  cathode  surface.  The  grid  was 
simulated  using  the  technique  of  para.  5.4.1  where  the  dc  potentials 
are  added  to  the  space  charge  potentials  to  obtain  the  total  potential  for 
each  iteration. 

In  order  to  test  the  effect  of  grids  on  the  noise  in  a  crossed-f ield 

(9) 

gun,  the  gridded  injected  beam  gun  used  by  Northrup  and  discussed  by  Fontana 
was  simulated  using  a  48  x  24  x  16  mesh  system.  Although  the  calculation  con¬ 
verged  in  approximately  50  iterations  an  error  in  the  dc  potentials  was  dis¬ 
covered  which  prevented  the  beam  from  entering  the  interaction  region. 
Unfortunately  time  did  not  permit  further  calculations  on  this  test  case. 
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Average  Beom  Temperature  (*K) 


Figure  6-5.  Bees  Temperature  In  the  x-z  Plane. 
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7.0  ACCURACY  AND  ERROR  CONSIDERATIONS  IN  COMPUTER  SIMULATIONS 

The  sources  of  error  due  Co  Che  computer  formulation  have  been 
Identified.  These  errors  can  be  divided  into  three  categories: 

(1)  Truncation  errors  due  to  the  finite  difference  approxi¬ 
mations  used  to  represent  the  Lorentz  Force  Equation 
(trajectory  equations)  and  Poisson's  Equation  (potential 
equation).  Such  errors  are  dependent  upon  the  magnitude 
of  the  mesh  intervals  in  both  time  and  space.  In  order 
to  minimize  these  errors »  the  finite  difference  equations 
have  been  carried  out  to  second  order.  The  convergence  of 
the  calculations  Indicates  that  the  finite  difference 
procedure  is  stable. 

(2)  Round-off  errors  due  to  the  Inexact  representation  of 
floating  point  numbers  in  the  computer.  The  error  de¬ 
pends  on  how  a  number  is  "rounded  off"  in  the  lowest  "bits" 
of  a  'word.'  The  three-dimensional  model  is  solved  on  an 
IBM  370/168  using  single  precision  floating  point  words 
which  have  an  accuracy  of  6  or  7  decimal  digits.  Clearly 
any  fluctuations  in  parameters  which  are  approximately  10**^ 
below  the  average  or  dc  value  cannot  be  calculated  accurately. 
Both  the  velocity  and  current  fluctuations  obtained  in  the 
model  are  approximately  two  orders  of  magnitude  larger  than 
the  round-off  values. 

(3)  Errors  Intrinsic  to  the  simulation  which  arise  from  the 
following  approximations: 

(a)  Representation  of  individual  electrons  by  a  point 
charge  (big  electron)  containing  many  electrons. 

This  approximation  results  In  an  exaggerated  collision 
frequency  for  the  electrons . 
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(b)  Representation  of  the  space  charge  assigned  to  a 
particular  mesh  point  by  a  uniform  cloud  centered 

on  a  mash  point  or  by  some  other  method.  This  error 
has  been  Investigated  by  Hockney  and  Birds  all. 

(c)  Representation  of  the  functions  by  finite  space  and 
time  steps  which  lead  to  a  lack  of  energy  conservation 
in  the  modal  and  thus  a  "heating"  of  the  electrons. 

This  problem  has  been  investigated  by  Hockney. 

The  magnitude  of  the  errors  caused  by  (a)«  (b) ,  and  (c) 
above  are  difficult  to  estimate.  Hockney ^^'^^as  investi¬ 
gated  various  techniques  which  reduce  these  errors.  Unfor¬ 
tunate  Ay,  the  techniques  suggested  by  Hockney  all  require 
Increased  computer  time  and  thus  Increase  the  cost  of  running 
the  program.  One  reasonably  simple  technique  called  the 
cloud-ln-cell  (C1C)  approach  has  proved  to  be  effective  in 
reducing  the  noise  due  to  (b)  and  (c)  above.  Unfortunately, 
time  did  not  permit  the  Implementation  of  the  cloud-in-cell 
approach. 

(4)  Errors  in  the  model  due  to  the  approximation  of  the  exact 
electrode  shapes  with  planar  electrodes. 

In  the  two-dimensional  model  sketched  in  Figure  7-1  the  space 
charge  contribution  to  the  electric  field  is  found  using  the  Idealized 
rectangular  geometry  shown  by  the  dotted  lines.  This  results  in  errors 
since  some  of  the  charges  are  not  imaged  in  an  electrode  at  the  proper 
distance  from  the  charge. 

An  estimate  of  this  error  can  be  found  by  solving  for  the  electric 
field  due  to  a  line  charge  imaged  in  a  perfect  conductor  as  sketched  in 
Figure  7-2.  For  this  case,  the  electric  field  is  given  by 


(7-1) 
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Figure  7-1  •  Exact  Electrode  Shapes-Long  Kino  Gun  Design 


Figure  7-2.  A  Line  Charge  Parallel  to  a 
Conducting  Plane. 
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where  q  is  the  line  charge  and  d  is  the  distance  between  the  line 
charge  and  the  conductor.  Note  that  the  second  tern  is  due  to  the 
conductor  and  reduces  the  field  below  that  of  a  "free"  line  charge. 


If  d  is  in  error  by  an  amount  Ad  the  above  equation  becones 


1 _ 1 

x-d— Ad  x+d+Ad 


(7-2) 


The  percent  error  due  to  the  incorrect  position  of  the  conductor 
is  given  by 


E  -E 
*1  *0 

Electric  Field  Error  ■  — — x  100Z 


(7-3) 


Using  Equation  (7-1)  and  (7-2)  in  Equation  (7-3)  we  obtain 


Electric  FUld  Error  -  *  1001 


where  r  ■  x-d,  the  distance  from  the  line  charge  to  the  field  point, 
is  kept  fixed  and  A  *  Ad/d.  This  equation  is  plotted  in  Figure  7-3 
as  a  function  of  r/d  for  various  values  of  A  as  parameters. 

tfe  note  from  Figure  7-3  that  as  long  as  the  field  point  is 
relatively  close  to  the  charge  (r/d  «  1),  the  error  in  the  value  of 
the  electric  field  is  small  even  for  large  errors  in  the  position  of  the 
conducting  plane  (A  •  1). 


For  the  gun  of  Figure  7-1  the  top  of  the  beam  is  approximately 
20  mesh  units  from  the  exact  electrode  and  28  mesh  units  from  the  approxi¬ 
mate  electrode.  This  gives  A  -  .4.  The  beam  thickness  is  approximately 
15  mesh  units.  Thus,  at  the  middle  of  the  beam  r/d  “  .33  and  the 
error  in  the  electric  field  at  the  center  of  the  beam  would  be  only  4.2Z. 
In  fact,  the  error  in  the  total  field  at  the  center  of  the  beam  would  be 
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chorqe  to  point  distonce _ 

r  charge  to  conducting  plane  distonce 

Figure  7-3.  Electric  Field  Error  Due  to  Incorrect 
Electrode  Position. 
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even  less  since  the  field  et  the  center  of  the  bean  would  be  dominated 
by  near  charges  rather  than  those  at  the  edge  of  the  beam. 

A  larger  error  would  occur  at  the  bottom  of  the  beam  where  at 
exit  the  beam  la  approximately  6  mesh  units  from  the  exact  electrode  and 
11  mesh  units  from  the  approximate  electrode.  This  gives  a  value  of 
A  ■  .833  and  at  the  center  of  the  beam  r/d  "  1.2.  From  Figure  7-3 
the  error  Is  approximately  20. 5Z. 

Again,  the  actual  error  would  be  smaller  due  to  the  shielding 
effects  of  the  nearer  charges. 

It  should  be  noted  that  In  the  critical  region  of  the  cathode 
surface,  the  Idealized  and  exact  gun  geometries  coincide  so  there  Is 
negligible  error  In  the  calculation  of  electric  field  near  the  cathode 
surface.  The  field  due  to  charges  Is  not  the  total  field  since  there 
Is  a  potential  applied  between  the  accelerator  and  cathode  electrodes. 
This  would  further  reduce  the  net  error  In  the  electric  fields. 


8.0 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  shows  that  it  Is  possible  to  aodel  three-dimensional 
crossed-fleld  beams  using  computer  simulations  and  to  obtain  useful  Infor¬ 
mation  about  the  beam  characteristics  from  the  models. 

Further  studies  of  the  three-dimensional  beams  with  grids  over  the 
cathode  would  be  useful  in  order  to  determine  the  effect  of  grids  on  the 
beam  noise  and  optics. 
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A.l 


Introduction 


This  survey  of  papers  on  noise  and  noise  mechanisms  in  crossed- 
field  electron  beams  is  not  Intended  to  be  comprehensive,  but  rather  is 
Intended  to  be  definitive  in  establishing  the  current  state  of  knowledge 
concerning  noise  mechanisms  in  crossed-field  beams. 

Three  noise  mechanisms  have  been  proposed  to  account  for  the 
noise  growth  in  crossed-field  beams:  diocotron  gain,  feedback  to  the 
potential  minimum  and  noise  enhancement  due  to  electron  travel  along 
magnetic  field  lines  (three-dimensional  effects).  Each  of  these 
mechanisms  is  discussed  below,  and  an  effort  is  made  to  establish 
current  knowledge  concerning  each  mechanism,  noise  measurements  re¬ 
lated  to  each  mechanism  and  possible  methods  of  noise  reduction. 

A.  2 _ Diocotron  Noise  Growth 

Below,  theoretical  calculations,  experimental  measurements 
and  possible  methods  of  reducing  diocotron  noise  growth  are  cited  from 
the  available  literature.  It  should  be  noted  that  accurate  and  de¬ 
finitive  quantitative  measurements  are  difficult  to  make,  which  explains 
in  part  the  lack  of  experimental  measurements.  It  should  also  be  noted 
that  the  theoretical  calculations  cited  are  only  approximate  since  the 
complexity  of  the  problem  makes  it  very  difficult  to  obtain  closed  form 
solutions  even  in  the  small  signal  regime  where  most  of  the  analyses 
apply.  One  other  problem  with  the  theoretical  and  experimental  results 
which  should  be  cited  is  the  almost  total  lack  of  noise  data  on  high 
current  density  "long"  guns  since  most  of  the  experimental  laboratory 
work  was  done  on  low  current  density  beams. 


A. 2.1  Theoretical  Calculations  and  Experimental  Results 

In  an  analysis  of  vave  propagation  In  slipping  streams  of 
electrons,  MacFarlane  and  Bay  (1949)  found  that  fluctuations  in  a 
stream  of  electrons  will  grow  even  without  the  presence  of  a  alow 
vave  circuit.  This  effect  in  crossed-field  beams  has  come  to  be  known 
as  the  "diocotron"  affect.  Gould  (1957)  developed  an  equation  for 
the  growth  rate  of  signals  due  to  this  effect  between  parallel  con¬ 
ducting  planes  for  both  centered  and  non-centered  beams.  Later,  Sasaki 
and  Van  Duzer  (1966)  used  coupled  mode  theory  to  obtain  essentially 
the  same  result  as  Gould.  The  results  of  Gould  and  Sasaki  and  Van  Duzer 
were  limited  In  that  only  the  three  space  charge  waves  were  used  In  the 
analysis,  whereas.  In  general,  there  are  five  possible  waves  which  can 
propagate  on  a  crossed-field  beam.  Finally,  Mantena  (1968)  considered 
aU  five  waves  and  developed  a  general  transfunction  matrix  which  gives 
the  diocotron  growth  for  a  thin,  non-accelerating  beam  In  a  drift  region 
between  parallel  conducting  boundaries.  Unfortunately  none  of  the  above 
results  can  be  applied  to  the  gun  region  since  the  beam  Is  accelerating 
and  the  electrodes  In  the  gun  region  are  not.  In  general,  parallel. 
However,  numerous  experiments  have  confirmed  that  the  growth  of  waves 
due  to  the  diocotron  effect  does  take  place  and  the  theoretical  re¬ 
sults  for  thin,  non-accelerating  beams  in  the  small  signal  regime 
agree  closely  with  experiment. 

Van  Duzer  (1961)  was  the  first  one  to  try  to  calculate  the 
growth  of  fluctuations  In  the  gun  region.  Using  modified  Utiwellyn- 
Peterson  equations  he  was  able  to  approximately  calculate  the  growth 


rat*  of  velocity,  current,  and  position  fluctuations  along  the  beam. 
Later,  Van  Duzer  (1963)  used  these  results  to  try  to  calculate  the 
noise  figure  of  a  crossed  field  amplifier  Including  a  gun  region,  a 
drift  region  and  an  amplifying  region.  His  calculated  valuea  were 
very  much  higher  than  those  measured  on  an  amplifier  operating  at  a 
frequency  of  1874  MHz.  Using  an  eigenvalue  perturbation  technique 
together  with  Monte  Carlo  calculations  Wadhva  (1964)  concluded  that 
the  growth  of  fluctuations  in  the  beam  is  mainly  due  to  diocotron 
growth  during  the  latter  part  of  the  beam  trajectory  in  the  gun  region. 
He  also  cited  some  experimental  work  which  seems  to  confirm  the  above 
conclusion.  Using  the  results  of  his  calculations  in  1964  for  the  gun 
region  Wadhwa  (1968)  calculated  the  noise  figure  of  an  injected  beam 
CFA  Including  both  a  three  and  five  wave  analysis.  A  comparison  of 
his  results  with  those  of  Mantena  and  Van  Duzer  (1963)  shows  good 
correlation  between  experiment  and  theory. 

A. 2. 2  Reduction  of  Diocotron  Noise  Growth 

In  order  to  reduce  the  growth  of  fluctuations  due  to  the 
diocotron  effect,  various  schemes  have  been  proposed.  These  tech¬ 
niques  are  reviewed  below  together  with  an  attempt  to  evaluate  the 
uaefulness  of  each  approach. 

A.2.2.1  Modification  of  beam  and  geometry  design  parameters 

By  examining  the  gain  equation  for  the  diocotron  growth, 

Wadhva  (1964)  proposed  several  design  conditions  which  would  reduce 
the  diocotron  growth.  A  low  current  density,  high  velocity,  high 


■agnatic  flald  baam  would  balp  raduca  tba  diocotron  growth.  Unfortunately , 
for  practically  all  of  the  application  now  envisioned  for  CFA's  rela¬ 
tively  low  velocity,  high  currant  density  beans  a^a  desired  with  asgnetlc 
fields  as  snail  as  possible  in  order  to  reduce  the  weight  of  the  deviue. 

Two  other  suggestions  of  Vadhwa  do  have  nerlt.  First,  he 
proposed  eliminating  the  length  of  the  transition  region  between  the 
gun  and  the  interaction  region.  This  has  been  done  in  all  current  CFA 
designs  and  although  no  careful  quantitative  neasurements  have  been 
aade  to  verify  that  the  noise  is  reduced,  qualitative  results  would 
indicate  that  it  is. 

Second,  he  proposed  that  the  beam  be  positioned  asymmetrically 
between  the  focusing  electrodes  in  the  gun  and  transition  regions.  This 
is  a  design  parsmeter  which  is  somewhat  under  the  control  of  the  designer. 

A.2.2.2  Noise  transformer 

Wadhva  and  Rowe  (1963)  proposed  that  a  noise  transformer  some¬ 
what  similar  to  that  used  in  O-type  tubes  night  be  useful  in  reducing 
the  noise  output  of  a  crossed  field  gun.  Wadhva  and  Van  Duzer  (1965,  1968) 
used  a  noise  transformer  in  an  experimental  S-band  CFA  and  achieved  a 
3.5  db  noise  figure.  They  attribute  the  low  noise  figure  obtained  to 
be  due,  in  part,  to  the  noise  transformer.  Sldhu  and  Wadhva  (1970) 
proposed  a  noise  transformer  composed  of  two  parallel  plate  regions 
and  suggest  that  such  a  device  between  the  gun  and  the  interaction 
apace  would  reduce  the  noise  as  well  as  allowing  wide  dynamic  range 
operation  of  the  gun. 


Van  Duzer  and  Barrla  (1964)  found  an  instability  in  trajectory  calcu¬ 
lation  they  ware  attempting  to  make  on  a  long  Kino  Gun  and  alao 
propoaed  a  nolae  mechanlem  caused  by  feedback  to  the  potential 
to  account  for  thla  apace  charge  Instability.  Davies  and  Polter  (1965) 
made  noise  measurements  on  a  long  Kino  Gun  and  were  able  to  Identify 
two  types  of  nolse-dlocotron  noise  and  noise  which  could  be  associated 
with  feedback  to  the  potential  minimum.  Hunn  and  Smol  (1964)  also 
found  that  the  noise  In  crossed-fleld  guns  Increases  with  the  length 
of  the  cathode. 

Bo  and  Van  Duzer  (1968)  proposed  a  feedback  model  for  the 
potential  minimum  region  of  a  crossed-fleld  gun  which  gave  at  least 
qualitative  agreement  with  the  experimental  results  cited  above.  They 
developed  a  stability  criteria  for  avoiding  excess  noise  due  to  feed¬ 
back  to  potential  minimum  of 

*k  <  -“c 

where  £.  Is  the  normalized  cathode  length  and  l  Is  the  normalized 
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length  of  a  cycloid  In  the  gun  region. 

It  seems  relatively  certain  that  the  mechanism  of  feedback 
to  the  potential  minimum  Is  causing  excess  noise  In  crossed-fleld  guns 
even  though  very  few  quantitative  measurements  have  been  made  and  the 
theory  due  to  Ho  and  Van  Duzer  is  only  approximate  and  limited  to  the 
small  signal  regime. 


A.2.2.3  Shielded  cathode 


Theoretical  work  by  Buneman,  Levy,  and  Line on  (1966)  Indicates 
that  a  partially  shielded  gun  (shielded  from  the  magnetic  field)  aay 
have  lower  dlocotron  growth  than  a  gun  where  the  nagnetlc  field  does 
not  lace  the  cathode.  No  experimental  work  has  been  reported  In  the 
area  of  shielded  cathodes. 

A.2.2.4  Non-equlpotentlal  cathodes 

PBtzl  (1966)  has  done  some  theoretical  work  on  non-equlpotentlal 
cathodes  which  Is  Interesting  from  a  theoretical.  If  not  practical,  point 
of  view.  His  calculations  Indicate  that  If  it  were  possible  to  build  a 
non-equlpotentlal  cathode,  a  more  stable  beam  would  result  with  decreased 
dlocotron  noise  growth.  No  experimental  work  has  been  attempted  In  the 
area  of  non-equlpotentlal  cathodes. 

A.  3 _ Peedback  to  the  Potential  Minimum 


A. 3.1  Experimental  Measurements  and  Theoretical  Calculations 

As  crossed-fleld  guns  were  made  "longer"  In  order  to  give 
higher  current  density  beams,  certain  anamolous  effects  were  observed. 
Epsztein  (1957)  noticed  that  large  values  of  sole  current  occurred  at 
certain  values  of  the  magnetic  field  In  "long"  guns.  Amaud  (1964)  was 
able  to  relate  these  results  to  the  ratio  of  the  height  of  the  beam 
above  the  cathode  to  the  length  of  the  cathode  and  proposed  that 
feedback  to  the  potential  m-tn-tmum  could  be  the  cause  of  the  excess  noise. 
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A. 3. 2  Reduction  of  Kolse  Growth  Due  to  Feedback  to  the  Potential  Minimum 

Since  feedback  to  the  potential  minimum  le  a  space  charge  effect 
any  scheme  or  method  vhich  reduced  the  space  charge  In  the  gun  region 
should  also  reduce  the  nolee  due  to  feedback.  Significant  reductions 
In  noise  were  observed  by  Arnaud  and  Doehler  (1961)  and  Amaud,  Diamond, 
and  Epszteln  (1962)  when  a  grid  was  added  above  the  gun  with  the  grid 
wires  perpendicular  to  the  direction  of  the  magnetic  field  and  the  grid 
close  to  the  cathode  surface.  One  possible  explanation  of  the  reduction 
In  noise  due  to  the  grid  is  the  reduction  of  space  charge  in  the  potential 
minimum  (fields  from  the  beam  terminating  on  the  grid  and  thus  reducing 
the  space  charge  fields)  and  thus  a  reduction  In  the  feedback  noise 
mechanism.  Such  a  conclusion  must  be  considered  tentative  since  no 
experimental  or  theoretical  work  has  baen  done  to  determine  why  or  how 
a  grid  reduces  the  noise. 

A.  4 _ Three-Dimensional  Effects 

All  of  the  theoretical  analyses  which  have  been  made  of  crossed- 
fleld  beams  have  been  two  dimensional.  Variations  and  fluctuations  of 
the  beam  In  the  magnetic  field  direction  have  been  Ignored.  Gandhi  (1965) 
made  some  calculations  which  Indicated  that  random  emission  velocities 
In  the  direction  of  the  magnetic  field  may  contribute  substantially  to 
the  noise.  Subsequently,  Slsodla  and  Gandhi  (1966),  Sisodla,  Gandhi,  and 
Kadhva  (1967)  and  Sisodla  and  Wadhwa  (1968)  found  that  large  reductions 
In  the  noise  output  of  a  crosaed-fleld  gun  could  be  achieved  by  either 
tilting  the  cathode  slightly  In  the  magnetic  field  direction  or  by 
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•lightly  perturbing  the  aagnetic  field  in  the  gun  region  to  give 
—11  aegnetlc  field  components  perpendicular  to  the  gun.  The  noise 
reduction  is  believed  to  be  due  to  smoothing  out  of  the  random  velo¬ 
cities  in  the  magnetic  field  direction,  although  this  conclusion  must 
be  considered  tentative  until  further  three-dimensional  theoretical 
•nd/or  experimental  work  la  completed. 
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ABSTRACT 

la  this  paper  a  self-consistent  computer 
aldad  CFA  gun  analysis  program  Is  discussed 
which  can  ba  ustd  to  solve  for  the  character- 
la tics  of  "long"  croasad  flald  guns  with 
arbitrary  alactroda  a hap as  and  including  tha  " 
af facta  of  apaca  charge.  Iba  raaulta  of  ^ 

analysing  a  "long"  Kino  gun  are  praaantad 
which  show  excellent  agzaanant  batwaan  tha  ‘  n*z 
thaoratlcal  gun  par  ana  tars  and  thoaa  obtainad 
by  tha  computer  analysis.  As  an  example  of 
tha  design  procedure,  tha  raaulta  of  deal going 
a  CFA  gun  and  transition  region  with  tha  accel¬ 
erator  at  anode  potential  are  glean.  Tha 
utility  of  tha  modal  la  further  11 lust rated 
by  tha  results  of  soma  prellalnsry  noise  cal¬ 
culations  made  by  simulating  tha  noise  charac¬ 
teristics  of  tha  beam  at  tha  cathode  surface. 


INTRODUCTION 

In  order  to  build  officiant,  high  power  in¬ 
jected  beam  crossed  flald  a^llflers  It  Is  nec¬ 
essary  to  ba  able  to  design  convergent  croasad 
field  guns  which  produce  relatively  thin  laminar 
beams.  Tha  deslpi  and  analysis  of  guns  for 
practical  CFA' a  is  complicated  due  to  the  cy¬ 
cloidal  notion  of  tha  electrons  which  leads  to 
trajectory  crossings  sad  thus  vary  complicated 
flow. 

In  tha  past  CFA  guns  have  bean  designed 
using  eapirlcal  techniques  or  the  design  methods 
suggested  by  Klno(l)  with  potential  plotting 
used  to  match  aqui potential  lines  In  tha  trans¬ 
ition  region  between  tha  gun  and  the  interaction 
region. 

A  self-consistent  CFA  gun  analysis  pro¬ 
cedure  Is  discusssd  below  which  can  ba  used  to 
solve  for  the  eba recta rl a ties  of  "lerg"  crossed 
field  guns  with  arbitrary  electrode  shapes • 

Tbs  ef facts  of  space  charge  ere  Included  in  the 
analysis.  In  addition,  the  analysis  procedure 
discussed  below  shows  prowls#  of  bsing  able  to 
be  used  to  study  noise  mechanisms  end  nolee 
growth  In  the  gun  and  transition  ragion. 


TWO  DIMENSIONAL  LAGRAMGXAR  MODEL 

In  this  approach  a  large  number  (up  to  6000) 
rode  of  eharge  are  used  to  simulate  the  beam  In 
the  gun  end  transition  regions  with  rods  being 
released  from  the  cathode  consistent  with  the 
emission  characteristics  of  the  cathode  surface. 

The  approach  la  discussed  below: 

(1)  First,  the  gum  end  transition  region  Is  * 
divided  into  a  rectangular  mesh. 

(2)  Next,  the  potentials  in  the  absence  of 
apaca  charga  are  calculated  at  each  mesh 
point  using  a  computer  code  which  solves 
LaPlace' s  aquation  for  arbitrary  alac¬ 
troda  shapes  and  potentials.  These 
potentials  are  then  stored  for  uee 
during  the  rest  of  the  calculation. 

(3)  A  relatively  large  number  of  rods  of 
charge  are  released  from  the  cathode 
with  position,  dsnsity,  and  velocity 
consistent  with  the  emission  charac¬ 
teristics  of  the  cathode.  The  haem 
buildup  is  followed  in  time  with  elec¬ 
trons  being  emitted  from  the  cathode  at 
equal  small  In  ere  man  ts  in  tins.  Tha 
beam  la  allowed  to  build  up  to  several 
thousand  rods  of  charge  with  "new"  rods 
of  charge  added  at  tha  cathode  In  each 
time  step  and  "old"  rods  being  removed 
aa  they  exit  tha  right-hand  boundary  or 
are  collected  on  the  electrodes. 

(A)  The  motion  of  the  rods  of  charge  is 

governed  by  the  Lo rents  Force  Equation 
where  the  electric  field  Is  found  from 
the  potentials  on  the  mesh  points. 

During  each  Iteration,  each  rod  of 
charge  la  allowed  to  move  for  a  small 
Increment  In  time  with  Its  motion 
governsd  by  the  Lorents  Force  Equation. 
Vary  accurate  finite  difference  aqua¬ 
tions  are  used  which  have  previously 
bean  developed  and  used  successfully 
to  solve  the  Lorants  Force  lquatlon(2). 
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(5)  Simulating  tba  btn  by  a  larga  number  of 
coda  of  charge  has  sot  boon  mid  1a  tba 
Hit  aa  an  approach  to  tba  daalgn  of 
•  croaaad  field  guna  dua  to  tba  relatively 
long  coop u tar  tlnaa  required  to  aolva 
Poisson's  equation  by  eonvaatlonal  tach- 
-  miquee.  However,  a  tachnlqua  dua  to 
L  Hockney (3)  baa  boon  auccaaafully  uaad 
to  solve  Poisson' a  Equation  In  laaa  than 
ana  aacond  on  aa  ISM  360/166  (cohered 
with  approximately  ona  nlnuta  for  tba 
oaual  Iterative  relaxation  procedures). 
Thla  approach  is  dlacuaaad  In  aoaa  da- 
tall  In  a  papar  by  Yut  Kooyara  and 
Sunaaan(2)  and  was  uaad  auccaaafully 
vltb  up  to  four  thousand  rods  of  eharga 
la  analyzing  for  tha  first  tins  tha  die- 
trlbutad  aalssloo  CPA.  Ones  tha  potan- 
tlala  hava  baan  calc ula tad  including 
apaca  eharga  by  tba  Hockney  tachnlqua 
for  each  Iteration*  "new"  rods  of  eharga 
are  released  from  tha  cathode  and  tha 
procedure  dlacuaaad  above  la  repeated 
until  the  baan  reaches  a  "steady  state." 
Steady  state  or  convergence  ia  obtained 
vhan  tha  nuaber  of  rods  of  charge  exit¬ 
ing  tba  gun  region  averaged  over  ona 
cyclotron  period  approaches  a  constant. 


GUM  STUDIES  USING  THE  MODEL 
Lons  Kino  Gun  Analysis 

In  order  to  teat  tha  accuracy  and  validity 
of  the  nodal  a  long  lino  gun  vaa  designed  with 
the  following  paranatars: 

Magna tic  Flald  -  .230  w/n2 

Cathode  Length  -  .130  In. 

Normalized  Length  -  116.4  Kino  units 

Accelerator  Potential  -  3000  volts 
Cathode  Currant  Danelty  -  14  amp/cm^ 

Total  Seen  Current  -  3.34  aap. 

Going  a  time  atap  of  1/20  of  a  cyclotron 
period,  the  been  buildup  charecterletlco  were 
calculated  for  110  Iterations.  Aa  rods  of  charge 
vara  injected  at  the  cathode  aurfaca  tha  bean 
current  aa  aeaoured  by  tha  rode  of  charge  exiting 
Che  gtm  region  rapidly  built  up  to  an  average 
value  of  3.24  eaps  which  differs  from  tha  Kino 
predicted  value  by  lose  than  32.  The  been  thick¬ 
ness  and  current  density  of  tha  exiting  baan  are 
also  in  close  agreeatnt  with  tha  theoretical 
values  given  by  tha  Kino  theory. 

Gim  Paslm  with  Accelerator  at  Cathode  Potential 

The  computer  modal  becomes  moat  useful  aa  a 
design  tool  for  guna  whan  uaad  to  design  both  tha 
gun  and  transition  regions,  for  a  particular 
duel-mode  CPA  a  gun  vaa  required  to  operate  at  a 
magnetic  flald  of  .3  v/m*  and  produce  a  current 
of  5  opt. 


) 


As  a  starting  point,  a  Kino  gun  daalgn  was 
made  with  the  following  characteristics: 

2 

Current  Density  -  14  aaps/ce 
Cathode  Length  -  .130  Inches 

Kino  Calculated  Current  -  3.7  amps 
Anode  Potential  -  6.2  kv. 

turning  tha  accelerator  at  full  anode  poten¬ 
tial  represents  aa  Innovation  made  possible  by 
recognising  that  In  a  grldded  gun  (this  design  w#~ 
aide  for  the  eventual  inclusion  of  a  grid) , 
grid  potential  can  ha  uaad  to  vary  tha  catbo 
currant.  Thus,  a  saps retd  accelerator  is  not 
oecesaary.  Eliminating  tha  accelerator  resu  <n 
three  limedlata  advantages: 

(1)  Seduction  of  tha  nu^er  of  raqulrad 
electrode  voltages  by  one. 

(2)  Elimination  of  tha  gap  between  the 
accelerator  end  the  anode,  thus  reducing 
the  length  required  for  the  transition 
region  which  in  turn  reduces  diocotron 
growth  and  gun  noise. 

(3)  Elimination  of  tha  potential  and  flald 
matching  problem  dua  to  the  accelerator 
and  anode  gap. 

Although  tha  cathode  current  fluctuates,  the 
average  value  given  by  the  coeputer  model  is  about 
5.0  amps  which  ia  only  slightly  lass  than  that 
predicted  by  tha  Kino  theory.  However,  it  oust  be 
kept  In  mind  that  the  accelerator  electrode  is  not 
the  exact  Kino  shape  due  to  the  necessary  com¬ 
promises  made  to  natch  up  the  anode  and  the 
accelerator  end  thus  it  la  not  expected  that  tha 
exact  Kino  current  would  be  obtained. 

The  beam  thickness  varies  with  position,  but 
sn  overage  value  of  about  10  ells  la  obtained  for 
tha  exiting  beam.  Thla  design  was  reached  only 
after  several  sole,  been  focusing  electrode,  and 
anode  shapes  were  analysed.  Only  with  e  coeputer 
model  which  la  capable  of  modeling  the  complete 
gun  region  beck  to  the  cathode,  can  designs  such 
as  this  ha  made  with  any  confidence. 

Noise  Simulation  in  the  Crossed  Field  Seam 

In  order  to  study  the  noise  properties  of  the 
beam,  it  le  necessary  to  simulate  the  emission 
properties  of  the  hot  cathode  surface  as  accurately 
aa  possible.  For  the  computer  model  this  mans 
that  tha  following  pa  ram  tars  must  be  determined 
at  tha  cathode  aurfaca: 

(1)  Tha  number  of  electrons  emitted. 

(2>  The  emission  velocities  of  tba  electrons . 

(3)  The  dm  of  emission  of  the  electrons. 
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Tha  Mtk»4  for  calculating  ucb  of  tht  abovi  thrnn 
yantlriii  is  i«mistd  tolov: 

(X)  Emission 

'  Xt  can  bo  shown  that  tho  probability  of 
nsactly  8  oloctroaa  balag  omitted  fron 
o  given  cathode  apot  of  length  Ay 
daring  the  tine  interval  At  is  given  by 


f(S)  - 

X  ia  the  average  number  of  elec- 
trona  emitted  in  each  At.  This  ia 
tecopilsed  to  be  the  Poisson  Distri¬ 
bution. 

the  probability  that  S  or  fever  elec¬ 
trode  are  emitted  in  the  tine  interval 
At  la  given  by  the  cumulative  distri¬ 
bution  function 


i  1 


»(s)  -  l  f(S) 

s*o 


0  <  KS)  (  l 

To  determine  the  nuaber  of  olectrona 
omitted  in  e  tine  interval.  At,  a 
uniformly  distributed  random  number 
between  0  and  1  la  generated.  The 
number  of  electrons  emitted  la  S  if 
• 

F(S-l)  <I(  F(S) 

When  ft  a  F(0)  no  electrons  ere  emitted. 

(X)  Kmieelon  Velocities 

A  half -Hews  Ilian  distribution  la  as- 
timid  for  the  x  component  of  emission 
velocity  (the  component  perpendicular 
to  the  cathode).  The  x  velocity  is 
given  by 


where  la  e  random  nuaber  uniformly 
distributed  between  0  end  1,  T  la  the 
cathode  ta^erature,  k  is  Aoltrmenn'a 
constant,  and  m  ia  the  electron  mass. 

A  fnll-Mawellian  distribution  is  as- 
aiaed  for  alectron  vnlocltlae  in  tha  y 
direction  (perellal  to  the  cathode). 
The  relation  giving  this  velocity  ia 


T  -  (*?)'• rt-‘<V 


where  ft  is  a  random  nuaber  wlformly 
distributed  between  0  end  1. 

O)  Imleelon  Times 

In  order  to  avoid  calculating  tha  posi¬ 
tions  of  tha  alectrona  above  the  cath¬ 
ode  for  random  emission  tines,  it  is 
assumed  that  sll*$  of  the  electrons  in 
a  given  mesh  region  are  emitted  at  the 
seme  time  but  with  a  random  distribu¬ 
tion  over  the  arse.  Tha  amission 
points  are  then  given  by 

*  -  *0*^? 

where  ft  ia  a  random  nuaber  between  0 
end  1  and  the  sign  of  ft  is  also 
selected  randomly.  This,  to  an  ob¬ 
server  looking  at  one  point  on  the 
cathode  surface  the  electrons  will 
appear  to  come  off  that  apot  at  random 
times.  Provided  the  time  interval  At 
is  chosen  to  be  relatively  small,  the 
degree  of  approximation  Introduced 
should  be  email. 

Preliminary  calculations  of  beam  "noisiness" 
as  measured  by  velocity  variance  end  nolee  power 
spectral'  density  calculations  have  led  to  the 
following  preliminary  observations : 

(1)  Operation  of  a  CPA  gun  In  the  te^erature 
limited  regime  leads  to  a  marked  Increase 
ia  beam  nolee.  Observation  of  the  beam 
shape  end  trajectories  during  teaperature 
limited  operation  shows  s  deterioration 
of  the  gun  optica  which  may  explain  the 
Increased  nolee. 

(2)  The  nolee  grows  approximately  exponen¬ 
tially  with  distance  ia  the  direction  of 
beam  flow.  This  result  is  not  imexpected 
•lace  amplification  of  the  noise  due  to 
the  dlocotron  effect  would  be  expected 

to  be  approximately  exponential. 

(3)  The  upper  edge  of  the  beam  (nearest  the 
anode)  is  noisier  then  the  lower  edge  of 
the  beam. 
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Computer-Aided  Design  of  Electron  Guns 
for  Injected- Beam  Crossed-Fieid  Amplifiers 

ELDEN  K.  SHAW  and  G.  P.  KOOYERS 

Abstract -k  two-dimensional  set f- consistent  computer-aided  CFA 
gun  analysis  program  is  discussed  which  can  be  used  to  solve  for  the 
characteristics  of  “long”  crossed -field  guns  with  arbitrary  electrode 
shapes  and  including  the  effects  of  space  charge.  The  results  of  analyz¬ 
ing  a  “long”  Kino  gun  tfiow  excellent  agreement  between  the  theoretical 
guft  parameters  and  those  obtained  by  the  computer  analysis.  The 
model  shows  promise  of  also  being  useful  in  the  analysis  of  noise  in 
CFA  guns. 

I.  Two-Dimensional  Lagrangian  Model 

In  this  approach  a  large  number  (up  to  6000)  of  rods  of 
charge  are  used  to  simulate  the  beam  in  the  gun  and  transition 
regions  with  rods  being  released  from  the  cathode  consistent 
with  the  emission  characteristics  of  the  cathode  surface. 

The  approach  is  summarized  below: 

1)  First,  the  gun  and  transition  region  is  divided  into  a 
rectangular  mesh. 

2)  Next,  the  potentials  in  the  absence  of  space  charge  are 
calculated  at  each  mesh  point  using  a  computer  code  which 
solves  Laplace’s  equation  for  arbitrary  electrode  shapes  and 
potentials  These  potentials  are  then  stored  for  use  during 
the  rest  of  the  calculation. 

3)  A  number  of  rods  of  charge  are  released  from  the 
cathode  with  position,  density,  and  velocity  consistent  with 
the  emission  characteristics  of  the  cathode.  The  beam  buildup 
is  followed  in  time  with  electrons  being  emitted  from  the 
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Fig.  1 .  Exact  electrode  shapes  for  “long”  Kino  gun  design. 


Fig.  2.  Computer- calculi  ted  beam  current  exiting  the  gun  region  as  a  function  of  time. 


ACCELERATOR 


|  CATHODE  I 

Fijg  3.  Beam  density  profile  for  the  “long”  Kino  gun.  Notes:  l)  There  art  2355  rods  of  charge  simulating  the  beam.  2)  152  rods  of  charge  are  in 
jected  at  each  time  step;  only  an  average  of  21  exit  the  gun  region.  3)  The  current  exiting  to  the  nght  differs  from  the  Kino  calculated  value  by 
only  5  percent.  4)  Random  electron  thermal  velocities  are  included  in  this  calculation.  5)  Electron  shapes  are  from  Kino  gun  theory. 


cathode  at  equal  small  increments  in  time.  The  beam  is  al¬ 
lowed  to  build  up  to  several  thousand  rods  of  charge  with 
“new"  rods  of  charge  added  at  the  cathode  in  each  time  step 
and  "old"  rods  being  removed  as  they  exit  the  right-hand 
boundary  or  are  collected  on  the  electrodes. 


4)  The  motion  of  the  rods  of  charge  is  governed  by  the 
Lorentz  force  equation  where  the  electric  field  is  found  from 
the  potentials  on  the  mesh  points.  During  each  iteration, 
each  rod  of  charge  is  allowed  to  move  for  a  small  increment  in 
time  with  its  motion  governed  by  the  Lorentz  force  equation. 
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Fig.  4.  Calculated  current  density  at  gun  exit  compared  to  theoretical 
beam  density  profile. 
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5)  Simulating  the  beam  by  a  large  number  of  rods  of  charge 
has  not  been  used  in  the  past  as  an  approach  to  the  design  of 
crossed-field  guns  due  to  the  relatively  long  computer  times 
required  to  solve  Poisson’s  equation  by  conventional  tech¬ 
niques.  However,  a  technique  due  to  Hockney  (11  -(3]  has 
been  successfully  used  to  solve  Poisson’s  equation  in  less  than 
1  s  on  an  IBM  370/168  (compared  with  approximately  1  min 
for  the  usual  iterative  relaxation  procedures).  This  approach  is 
discussed  in  some  detail  in  a  paper  by  Yu,  Kooyers,  and 
Buneman  (4|  and  was  used  successfully  in  analyzing  for  the 
first  time  the  distributed  emission  CFA.  Once  the  potentials 
have  been  calculated  including  space  charge  by  the  Hockney 
technique  for  each  iteration,  ’’new”  rods  of  charge  are  released 
from  the  cathode  and  the  procedure  discussed  above  is  re¬ 
peated  until  the  beam  reaches  a  ’’steady  state.”  Steady  state 
or  convergence  is  obtained  when  the  number  of  rods  exiting 
the  gun  region  averaged  over  one  cyclotron  period  approaches 
a  constant. 


11.  analysis  of  a  Long  Kino  Gun 


In  order  to  determine  the  validity  of  the  computer  model, 
the  ’’long”  Kino  gun  (5)  shown  in  Fig.  1  was  analyzed.  The 
parameters  of  this  gun  are. 


total  beam  current 
cathode  length 
magnetic  field 
current  density 
accelerator  potential 
normalized  cathode  length 


3.34  A 
0.1 50  in 
0.2500  W/m2 
14  A/cm2 
5000  V 

1 16.4  Kino  units. 


Current  exiting  ‘he  gun  region  as  a  function  of  time  step  is 
shown  in  Fig.  2.  For  this  case  the  current  rises  rapidly  and 
then  fluctuates  about  a  mean  value  which  is  close  to  the  Kino 
predicted  value.  In  this  case,  the  cathode  is  operating  space- 
charge-limited  with  only  a  fraction  of  the  total  emitted  rods 
exiting  the  gun  region.  Shown  in  Fig.  3  is  the  exact  Kino 
beam  shape  compared  to  that  obtained  by  plotting  the  in¬ 
dividual  charge  positions  from  the  computer  program.  Again, 
.he  agreement  between  the  computer-calculated  results  and 
the  Kino  theory  is  good.  In  Fig  4  the  computer-calculated 
current  density  is  compared  to  the  ’’exact”  values  calculated 
using  the  Kino  equations. 
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APPENDIX  D 


SOLVING  POISSON'S  EQUATION 
IN  THREE  DIMENSIONS 


D.l 


Solving  Poisson’s  Equation  in  Three  Dimensions 
We  want  to  solve  Poisson's  Equation  in  the  region  sketched  in 
Figure  D-l  with  a  constant  potential,  on  the  z  boundaries. 

First,  the  potential,  is  subtracted  from  all  boundary 

potentials  since  this  results  in  simpler  boundary  conditions  (l.e.,  two 
boundaries  at  zero  potential  rather  than  one).  The  modified  boundary 
conditions  are  shown  in  Figure  D-l. 

A  coordinate  grid  system  is  Introduced  with  1  segments  in 
the  x- direction,  m  in  the  y-dlrectlon,  and  n  in  the  z-direction  as 
sketched  in  Figure  D-l. 

For  the  mesh  system  defined  above  the  boundary  conditions  are 
now  given  by 


*i,J,l  “  *i,j,n+l  "  0 


*l,j,k  “  "*C 


*1+1, j,k  “  *A  “  *C  “  *B 


(D-l) 


*1,0, k  "  *i,2,k 


*l,m,k  "  *l,mt2,k 


ii  -  o 

dn 


The  five  point  difference  approximation  to  Poisson's  Equation  is 


Ax7  *1-1, j,k  “  2*i,j,k  +  *i+l, j ,k  + 


IF  *i,J-l,k  "  2*i,J,k+  *i,J+l,k  + 


IF  *i»J,k-l  "  2*i,J,k  +  *i,J,k+l 


JLxlA 


(D-2) 


Figure  D-l.  Three-Dimensional  Geometry  with  Modified 
Bojudary  Conditions  for  the  Case  of 
Constant  Potential  on  z-Boundary. 


o>  |a> 


(D-5) 


Substituting  the  expression  for  $  into  the  three-dimensional 
five  point  difference  equation*  simplifying  and  equating  the  sine  com¬ 
ponents*  we  obtain  the  equation 

Pi,j,s+1  "  Ax7  *i-l,j,s+l  “  2*i,J,s+l  +  *i+l,J,s+l 

+  JL  —  -  2—  +  —  +  -2—  cog  JLS. 

Ay7  ♦i.J-^s+l  >i,J,s+l  ♦i.J+l.s+l  A*7  n 


*i,J.s+l 


(D-6) 


where  we  have  also  used  the  sine  expansion  for  p.  We  note  that  this 
is  just  the  finite  difference  approximation  to  the  Helmholtz  equation 
in  two  dimensions  for  each  value  of  s.  Thus*  in  order  to  solve  the 
three-dimensional  Poisson  Equation,  we  use  the  following  procedure: 


Cl) 


Calculate  p.  .  .  for  all  l,j  by  analyzing  p,  .  , 

If  J  t  It  J  » 

along  k. 


Pi.J,k*l 


t 


(D-7) 


with  the  p.  .  . 's  being  known. 

J  >  * 

This  will  require  applying  a  Fourier  Analysis  program 
Ct-l)x(m-l)xCn-l)  times.  To  accomplish  this  analysis,  an 
efficient  computer  program  called  FOUR67  Is  used. 


(2)  Knowing  the  p's  we  then  solve  the  two-dimensional 
Helmholtz  equation 


V2<fr  +  K2*  -  p 


for  each  of  the  n-1  values  of  K.  This  two-dimensional 
problem  Is  solved  using  a  computer  program  called  PWSCRT. 


In  order  to  solve  the  Helmholtz  equation  the  boundary  con* 
ditions  of  the  transformed  potentials  must  be  known.  We 
note  that  since 


*i,J.k  "  ~*C 

“d  ♦l+l.j.k  "  V*C 

then  ?  .  -  -  4>c  l  «ln 

i.J.k  C  %ml  n 


*°d  k  ■  J  ,ln 


(D-i 


CD—! 


The  boundary  conditions  at  J*1  and  j-ntfl  are 


|i  -  ° 

3y 

or  ♦i.O.k  “  ♦i,2»k 

*l*m,k  "  ^l,m+l,k 
which  give 

art-1, k  m  8^l,l,k 
3y  3y 


CD-: 


These  boundary  conditions  are  sketched  in  Figure  D-2  for 
the  k-th  row  In  the  z-dlrectlon. 


ve  then  use  the  inverse  Fourier  Transform 


APPENDIX  E 


FINITE  DIFFERENCE  EQUATIONS  OF 
MOTION  IN  THREE  DIMENSIONS 


r 


E.l  EQUATIONS  OF  MOTION 


In  order  to  find  the  motion  of  the  charges.  It  Is  necessary 
to  solve  the  equation 


4—  “  q(?  +  v  x  I)  . 


md£  ‘  q 


Assuming  both  electric  and  magnetic  field  components  In  the  x, 
y,  and  z  directions,  the  components  of  the  above  equation  become 


s.  -  J ax ♦*».-» y 


(E-l) 


4  (E  +  z  B  -xB) 
m  y  x  z 


4  (E  -  y  B  +  iB  ) 
m  s  J  z  y 


(E-2) 


(E-3) 


•Assuming  that  the  electric  and  magnetic  fields  are  constant  in 
the  Interval  of  Integration,  the  -ove  equations  can  be  Integrated  once 


to  obtain: 


i(t)  -  x(0)  -  4  [^t  +  (y(t)  - 
y(t)  -  y(0)  -  4  ^yt  +  (*(t)  - 


y(0))Bz  -  (z(t)  -  s(0))b2 

(E-4) 


s(0))Bx  -  (x(t)  -  x(0))B  I 

(E-5) 


z(t)  -  z(0)  -  4  Fe  t  -  (y(t)  -  y(0))B  +  (x(t)  -  x(0)  )B 

m  |_z  z  5 


(E-6) 


Using  the  Laplace  Transform  on  Eqs<  E— 4,  E— 5  and  E— 6,  we  obtain 
the  linear  equations 

E  i 

_  _  ,  X  ,  A  /  .  .  _  _  _  X  i  \  _L  _ 


-  oj  y  +  w  z  m  n  ~"T  +  —  (u>  zn  -  w  yn  +  irt)  +  xn 
czJ  cy  8  cy  0  cz^O  0  0 


(E-7) 
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nE*  [“ex*  t2 

ft-  (-1  +  Y  +  co«t) 

c  Lc  J 


T)E  b)  U  _2  (0 

— 5  —  (-1  +  ~  +  cost)  +  —Si  (t  -  slnt) 

u  ^  u  u  2  u 

c  c  c  c 


nE,  l«r_  m  .2  “ 

- -  (-1  +  ■=—  +  cost)  -  —Si  (t  -  sIht) 

w_  w  w  '  2  or 

c  lc  c  c  _ 


(E-10 


where  u  2 
c 


(l)  2  +  (1 i  2+u  2 

cx  cy  cz 


end  t  ■  o»  ,  . 

ct 


Equations  for  y  and  s  are  not  given  since  they  are  easily 
obtained  from  Equation  E-10  by  permuting  the  coordinates  as  given  below: 


*  ♦  y 

7  *►  * 
z  <+•  x 


In  the  finite  difference  procedures  used  in  the  computer  program 
the  x,  y»  and  s  positions  of  the  particles  are  found  at  t  +  At  using 
Equation  E-7  and  the  similar  equations  for  y  and  z.  Equations  E-4, 
E-5,  and  E-6  are  then  used  to  find  the  x,  y,  and  z  velocities  at 
t  +  At. 
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